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ABSTRACT 

We present direct observational evidence for star formation quenching in galaxy groups in the redshift 
range 0 < z < 2.5. We utilize a large sample of nearly 6000 groups, selected by fixed cumulative 
number density from three photometric catalogs, to follow the evolving quiescent fractions of central 
and satellite galaxies over roughly 11 Gyr. At z ~ 0, central galaxies in our sample range in stellar 
mass from Milky Way/M31 analogs (M*/M 0 = 6.5 x 10 10 ) to nearby massive ellipticals (M*/M 0 = 

1.5 x 10 11 ). Satellite galaxies in the same groups reach masses as low as twice that of the Large 
Magellanic Cloud (M*/M® = 6.5 x 10 9 ). Using statistical background subtraction, we measure the 
average rest-frame colors of galaxies in our groups and calculate the evolving quiescent fractions of 
centrals and satellites over seven redshift bins. Our analysis shows clear evidence for star formation 
quenching in group halos, with a different quenching onset for centrals and their satellite galaxies. 

Using halo mass estimates for our central galaxies, we find that star formation shuts off in centrals 
when typical halo masses reach between 10 12 and 1O 13 M 0 , consistent with predictions from the halo 
quenching model. In contrast, satellite galaxies in the same groups most likely undergo quenching by 
environmental processes, whose onset is delayed with respect to their central galaxy. Although star 
formation is suppressed in all galaxies over time, the processes that govern quenching are different for 
centrals and satellites. While mass plays an important role in determining the star formation activity 
of central galaxies, quenching in satellite galaxies is dominated by the environment in which they 
reside. 

Keywords: galaxies: evolution - galaxies: groups: general - galaxies: star formation 


1. INTRODUCTION 

Star formation activity in galaxies at low redshift 
varies strongly with local environment. Correlations be- 
tween local galaxy density and star formation tracers 
(e.g., star formation rates, stellar colors, morphology) 
have long been observed in the nearby universe (e.g., 
Oemler 1974; Dressier 1980; Kauffmann et al. 2004; 
Hogg et al. 2004; Blanton et al. 2005a; Thomas et al. 
2005). Many physical mechanisms have been shown to 
be relevant for regulating star formation in dense envi- 
ronments by driving cold gas away from galaxies and by 
heating it up. The high fraction of non star forming 
galaxies in clusters and massive groups implies that such 
mechanisms have already acted to quench star formation 
in low redshift halos. 

Tremendous effort has been devoted in recent years to 
analyses of galaxy groups and lower mass environments. 
While such studies find that the properties of galaxies in 
groups are indeed different than those of field galaxies 
(e.g., Zabludoff & Mulchaey 1998; Balogh et al. 2004, 
2009; Wilman et al. 2005; Weinmann et al. 2006; Cooper 
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et al. 2007; Tran et al. 2009; Patel et al. 2009, 2011; 
Bolzonella et al. 2010; Peng et al. 2010; Kovac et al. 
2010; McGee et al. 2011; Quadri et al. 2012; Muzzin et 
al. 2012; Geha et al. 2012; Rasmussen et al. 2012; Hou 
et al. 2013), it appears that environmental processes do 
not influence all galaxies equally (e.g., Yang et al. 2005; 
van den Bosch et al. 2008; Font et al. 2008; Skibba & 
Sheth 2009; Skibba 2009; Guo et al. 2009; Peng et al. 
2012; Wetzel et al. 2012, 2013; Knobel et al. 2013; Woo 
et al. 2013; Carollo et al. 2013). 

Central galaxies, defined as the most massive galax- 
ies in their group or cluster, are affected by processes 
whose star formation quenching efficiency is proportional 
to the mass of the surrounding host halo. Such processes 
can act on the central galaxy itself (e.g., active galac- 
tic nuclei: Kauffmann et al. 2000; Springel et al. 2005; 
Hopkins et al. 2006; Fabian et al. 2012) or they can 
influence galaxies on halo scales (e.g., halo- quenching: 
Birnboim & Dekel 2003; Dekel & Birnboim 2006). In 
addition to being affected by the same processes, satel- 
lite galaxies in groups and clusters are also subject to 
gas and stellar stripping, strangulation and harassment, 
all of which contribute to a gradual shut off of star for- 
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mation (e.g., Gunn & Gott 1972; Larson et al. 1980; 
Farouki & Shapiro 1981; Byrd & Valtonen 1990; Moore 
et al. 1996; Balogh et al. 2000). Such interactions take 
place between satellite galaxies or between satellites and 
the host halo itself, and are sensitive to the varying local 
gas and galaxy density in the halo. 

Internal processes are likely also important drivers of 
star formation quenching in all galaxies. Recent studies 
find that structural properties of galaxies, such as central 
mass density and bulge-to-total mass ratio, are correlated 
with star formation activity at low and high redshift (e.g., 
Cheung et al. 2012; Fang et al. 2013; Barro et al. 2013). 
In addition, stellar velocity dispersion was found to be a 
good proxy for galaxy color and star formation activity, 
essentially linking the inner stellar body with the host 
dark matter halo (e.g., Franx et al. 2008; Wake et al. 
2012; Bezanson et al. 2012; Weinmann et al. 2013). Such 
correlations indicate that the local properties of galaxies 
pose essential conditions for star formation suppression. 

Here we perform a statistical analysis of galaxies in 
groups at 0 < z < 2.5. Since groups assemble and evolve 
later than massive clusters, quenching processes are still 
important even at relatively late times. We are able 
to explore the average properties of central and satellite 
galaxies over a large redshift range using only photomet- 
ric catalogs and statistical background subtraction. This 
approach is useful for deriving the evolution of galaxy 
properties without the need of spectra (e.g., Masjedi et 
al. 2006, 2008; Watson et al. 2012; Tal et al. 2012, 
2013; Wang & White 2012; Marmol-Queralto et al. 2012, 
2013). This may currently be the most efficient way to 
study galaxy group halos even at moderate redshift. 

Throughout the paper we adopt the following cosmo- 
logical parameters: = 0.3, = 0.7 and Hq = 70 

km s -1 Mpc -1 . 

2. PHOTOMETRIC CATALOGS AND SAMPLE 
SELECTION 

Galaxies for this study were selected at 0 < z < 2.5 
from three datasets, utilizing data from four surveys: Ul- 
traVISTA (McCracken et al. 2012), 3D-HST (Brammer 
et al. 2012), Cosmic Assembly Near-infrared Deep Extra- 
galactic Legacy Survey (CANDELS; Grogin et al. 2011), 
and Sloan Digital Sky Survey (SDSS; York et al. 2000). 

Low redshift galaxies were selected at 0.02 < z < 0.04 
from the NYU Value-Added Galaxy Catalog (NYU- 
VAG; Blanton et al. 2005b), which was derived using 
imaging and spectroscopic data from the seventh data 
release of SDSS (Abazajian et al. 2009) and imaging 
data from the Two Micron All Sky Survey (Skrutskie 
et al. 2006). Redshifts were determined spectroscopi- 
cally as part of SDSS for roughly 94% of the galaxies in 
that catalog while the rest of the galaxies were assigned 
the same redshift measurement as their nearest spectro- 
scopic neighbor. All galaxies in our low redshift sample 
are more massive than 5 x 10 9 M Q , the stellar mass com- 
pleteness limit at z = 0.04 (e.g., Wake et al. 2012). 

Galaxies in the redshift range 0.2 < z < 1.2 were 
identified from a public K s -selected catalog (Muzzin 
et al. 2013b) based on the first data release of Ul- 
traVISTA, an ongoing ultra deep near-infrared survey 
with the European Southern Observatory VISTA sur- 
vey telescope. The catalog covers a total area of 1.62 
deg 2 in the COSMOS field. It includes photometry in 
30 bands and provides excellent photometric redshifts 
(cr z /(l + z) = 0.013). Stellar mass completeness analysis 



Figure 1 . Completeness in the 3D-HST/CANDELS catalog as a 
function of galaxy stellar mass. Completeness was calculated by 
comparing object detection in CANDELS/deep with a re-combined 
subset of the exposures which reach the depth of CANDELS /wide 
(see Section 2.1 for details). Each of the curves represents this 
measurement in a different redshift bin. The catalogs are 90% 
complete down to 1.1 X 1O 9 M0 at 2 < 1.8 and down to 2.6 X 10 9 M© 
at 2 < 2.5. 

was performed for this data set by Muzzin et al. (2013a), 
who found a 95% completeness limit of 5 x 10 9 M Q at 
z < 1.2. 

Galaxies at 2 > 1.2 were observed as part of the 
complimentary 3D-HST and CANDELS surveys. These 
are near-infrared spectroscopic and imaging surveys with 
the Wide Field Camera 3 on board the Hubble Space 
Telescope. Images and grism spectra are available from 
the surveys in five fields (AEGIS, COSMOS, GOODS-S, 
GOODS-N and UDS), covering a total area of roughly 
0.2 deg 2 . Photometric catalogs, including excellent pho- 
tometric redshifts (cr z /(l + z) ~ 0.02) and stellar popu- 
lation synthesis model fits were derived by the 3D-HST 
team using images from both surveys in addition to a 
large set of ancillary data (R. Skelton 2014, in prepara- 
tion). We determine the stellar mass completeness limits 
of the 3D-HST/CANDELS data set in the following sec- 
tion. 

2.1. Completeness in 3D-HST /CANDELS 

We estimate the completeness in stellar mass of the 
3D-HST/CANDELS v3.0 catalog by comparing source 
detection in the wide and deep components of the CAN- 
DELS survey. We do so by creating a detection image 
over the CANDELS/deep area with the same exposure 
time and depth as the shallower CANDELS/wide survey. 
Although CANDELS/deep covers only about 17% of the 
total survey area, it reaches roughly twice the depth of 
CANDELS/wide (with approximately four times the ex- 
posure time), thus making it suitable for calculating the 
fraction of missed objects as a function of stellar mass. 

Following the same methodology that was used to cre- 
ate the science grade images, we reduced only 25% of 
the data that is available in the GOODS-S/deep field to 
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Figure 2. Sample matching across different redshift bins based 
on constant number density. The color curves represent the cu- 
mulative stellar mass functions from Marchesini et al. (2009). We 
calculated galaxy masses at four constant cumulative number den- 
sity values and derived mass-redshift relations which we fitted with 
second order polynomials (inset figure). For comparison, the square 
data points in the inset figure show the same measurement using 
the stellar mass functions of Muzzin et al 2013a. The selection 
at n = 4 x 10 -4 Mpc -3 is marked for reference (gray horizontal 
line and thick curve in the inset figure). We utilize the resulting 
mass-redshift relations to select central galaxies at 0 < z < 2.5. 


Table 1 

Central galaxy selection criteria at n = 4 x 10 — 4 Mpc -3 


Data set 

Redshift 

Mmed/M© 

Rv [kpc] 

Ncen 

SDSS 

0.02 < 2 < 0.04 

7.9xl0 10 

250 

251 


0.2 < z < 0.4 

6.9xl0 10 

220 

261 

UltraVISTA 

0.4 < 2 : < 0.6 

6.1xl0 10 

200 

460 


0.6 < z < 0.85 

5.2xl0 10 

180 

1080 


0.85 < z < 1.2 

4.1xl0 10 

150 

2281 


1.2 < 2 : < 1.8 

2.5xl0 10 

110 

707 

3D-HST 

1.8 < z < 2.5 

1.2xl0 10 

75 

915 


2.5 < 2 : < 3.5 

2.8xl0 9 

40 

3054 


2.2. Cumulative Number Density Matching 

Several recent studies have demonstrated that select- 
ing galaxy samples at a fixed cumulative number density 
effectively allows tracing the evolution of a given popu- 
lation of galaxies over different epochs (e.g., Gao et al. 
2004; van Dokkum et al. 2010; Papovich et al. 2011; 
Leja et al. 2013; Behroozi et al. 2013a). In this study 
we followed the same technique to match galaxy pop- 
ulations over the redshift range 0.02 < z < 3.5. We 
calculated the cumulative number density curves of the 
redshift samples from Marchesini et al. (2009) and mea- 
sured the implied stellar mass evolution at four num- 
ber density values: n = 1,2, 4, 6 x 10 -4 Mpc -3 (Figure 
2). For example, the resulting mass-redshift relation at 
n = 4 x 10 -4 Mpc -3 is: 

log(M*/M©) = -0.10 z 2 - 0.18z + 10.90 (2) 


match the depth of the wide component. From CAN- 
DELS/deep we selected observations from visits 1 and 
2 (as defined in Grogin et al. 2011) and combined 
them to create 2-orbit depth images in F125W and 
F160W. From 3D-HST we combined all F140W direct 
images in GOODS-S into a single orbit depth image. Fi- 
nally, we created a noise equalized detection image fol- 
lowing the same procedure that was used to produce 
the 3D-HST/CANDELS catalog. We multiplied each 
of the three science grade images (scp) by the square- 
root of their corresponding exposure-time weight-mask 
(©/ weight-) and averaged them according to the follow- 
ing equation: 


Detection image = 


Y (x/weightj x scij 
Y y weighty 


( 1 ) 


Source identification and matching was then performed 
on both depth level images. We ran SExtractor (Bertin 
& Arnouts 1996) on the newly combined “shallow” image 
as well as on the deep GOODS-S image using the same 
parameters that were utilized to create the photometric 
catalog. Figure 1 shows the fraction of detected objects 
in the shallow image compared to the deep image as a 
function of mass and redshift. Input redshift and stellar 
mass estimates were taken from the 3D-HST / C ANDELS 
catalog, for which detected galaxies from both images 
were matched to within 3 pixels, the typical size of the 
WFC3 point spread function. On that basis we esti- 
mate that the catalogs are 90% complete in stellar mass 
down to log (M s /Mq) = 9.04 at z < 1.8 and they reach 
log(M s /M 0 ) = 9.41 at z < 2.5. At z < 3.5 the catalogs 
are 75% complete down to log (M s /M©) = 9.14. 


Similar relations for the other number density samples 
are given in Appendix A. Toy model based estimates of 
the mass evolution at z > 1 imply a similar slope at high 
redshift (Dekel et al. 2013). 

We note that while the cumulative number density 
matching technique provides an estimate for the mass 
evolution of a given galaxy population over time, it re- 
lies on the assumption that the rank order of galaxies 
does not change dramatically over the studied redshift 
range. The presence of a significant major merger rate 
may change the mass ranking of galaxies and therefore 
render this assumption inadequate. In addition, a high 
major merging rate would affect this measurement by di- 
rectly changing the number of galaxies at a given mass. 

2.3. Central Galaxy Identification 

We divided all galaxies from the UltraVISTA and 3D- 
HST/C ANDELS catalogs into six redshift bins such that 
each bin spans roughly 1.5 to 2 Gyr. The time span of 
the SDSS redshift bin is shorter, as it was limited by 
the catalog mass completeness limit. At each redshift we 
selected central galaxy candidates in log M* bins of width 
0.3 dex with an evolving median mass according to the 
above calculated mass-redshift relation (e.g., Equation 
2). We considered a galaxy to be central if no other, more 
massive, galaxy was found within two projected virial 
radii of the host halo of the central candidate. Virial 
radius estimates for host halos around centrals at a given 
median mass and a given redshift were measured from 
catalogs extracted from the semi-analytic model of Guo 
et al. (2011). Median stellar mass, redshift limits and 
virial radius estimates for central galaxies in each of the 
bins at n = 4 x 10 -4 Mpc -3 are given in Table 1. At 
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Figure 3. Redshift evolution of the fraction of quiescent central 
galaxies. Each set of points represent a fixed cumulative number 
density selection of central galaxies and are fitted by a line. Stellar 
mass estimates at z = 0 are also noted for each value of n. Error 
bars are calculated as propagated y/N noise estimates. While the 
fraction of quiescent galaxies increases with time for all samples, 
lower number density samples (which correspond to earlier mass 
growth) preferentially get quenched at higher redshift. 

z = 0, central galaxies in our sample reach M*/M 0 = 
6.5 x 10 10 , roughly the mass of the Milky Way and M31. 
Detailed selection criteria for all samples are presented 
in Appendix A. 

2.4. Rest frame UVJ Colors 

Global star formation rates have been measured in dis- 
tant galaxies using numerous techniques, often relying on 
photometric estimates of rest frame colors. Williams et 
al. (2009) showed that galaxies out to z = 2 can be reli- 
ably separated into either “star forming” or “quiescent” 
based on a rest frame U — V vs. V — J {UVJ) color 
selection. Recently, Muzzin et al. (2013b) extended this 
measurement to z = 3.5 using redshift and rest frame 
color estimates from the UltraVISTA survey catalog. 

Following the same approach, we identified each galaxy 
in our sample as either star forming or quiescent based on 
the derived rest frame colors in the catalogs. We adopted 
the UVJ threshold values from Muzzin et al. (2013b) to 
separate galaxies at z > 0.2 in the UltraVISTA and 3D- 
HST/CANDELS catalogs. At lower redshift we found 
equivalent thresholds in the u — g vs. g — J color plane 
that we used for selecting galaxies from the NYU-VAG 
catalog. We note that the relatively poor u— band pho- 
tometry in SDSS implies that rest frame u — g color es- 
timates at low redshift are significantly more scattered 
than U — V estimates in UltraVISTA and 3D-HST (see 
Appendix B for details). 

3. HALO QUENCHING OF CENTRAL GALAXIES 

Star formation suppression in central galaxies is often 
linked to the mass of their surrounding host halos. In this 
“halo quenching” model, strong shocks in halos above a 
critical mass M cr [ t ~ 10 12 M Q heat up any infalling cold 
gas and thus prevent further star formation in the central 



Figure 4. Predicted and observed redshift evolution of the thresh- 
old mass for halo quenching. Halo mass proxy estimates are calcu- 
lated under Press Schechter formalism for each pair of stellar mass 
and number density values. The red points and lines represent the 
halo mass estimates and redshifts of each number density sample 
when the fraction of quiescent galaxies reaches 10%, 20%, 30% 
and 50%. Also plotted are theoretical predictions from the halo- 
quenching model, as calculated by Cattaneo et al. (2006, magenta 
line) and by Dekel & Birnboin (2006, green line). The observed 
typical halo mass at which central galaxies are quenched is be- 
tween 10 12 Mq and 10 13 Mq, broadly in agreement with predicted 
values. 

galaxy itself (e.g., Birnboim & Dekel 2003; Dekel & Birn- 
boim 2006). This model predicts that the value of M crit 
is mostly constant at low redshift and that at high red- 
shift significant cold gas inflows allow star formation even 
in halos more massive than M cr i t . Recently, Behroozi 
et al. (2013b) matched observed galaxies to simulated 
dark matter halos to show that star formation efficiency 
peaks in halos with a typical mass around 5 x 10 n M o 
and that it drops quickly in more massive halos. Here 
we follow an alternative approach and estimate the halo 
mass at which quenching becomes important by analyz- 
ing the fraction of quiescent central galaxies as a function 
of redshift and mass. 

3.1. Early quenching of central galaxies 

We use rest frame UVJ color estimates to identify 
quenched central galaxies in each redshift sample. Figure 
3 shows the redshift dependence of the fraction of quies- 
cent central galaxies in each of the four number density 
selected samples. The red circles and line represent the 
n = 4 x 10 -4 Mpc -3 sample while lighter and darker 
points and lines show the same measurement at higher 
and lower fixed number density values. Also noted in 
the figure are z = 0 stellar mass estimates of galaxies in 
each sample. The fraction of quiescent central galaxies at 
all number densities seem to be increasing at nearly the 
same rate at z < 2.5. The evolving quiescent fractions 
that are presented in Figure 3 are consistent with results 
from recent similar analyses at fixed cumulative number 
density (Patel et al. 2013; Lundgren et al. 2013). 

Figure 3 shows that central galaxy quenching occurs 
early. The quiescent fraction of the lowest number den- 
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sity sample (1.5 x 10 n M© at z = 0) reaches 30% already 
at z > 2.0. Lower mass central galaxies, drawn from the 
higher number density samples, also begin quenching at 
relatively early times and are more than 10% quenched 
by z = 2.0. 

This result agrees qualitatively with the halo quench- 
ing model which predicts that centrals cease forming 
stars promptly after their host halo reaches the thresh- 
old mass M cr i t . In the following Section, we estimate the 
evolution of this threshold mass from observations. 

3.2. Halo Mass and Quenching 

We examine the influence of halo mass on star for- 
mation suppression in central galaxies by estimating the 
typical host halo mass at four fixed values of quenched 
fractions (f q = 10%, 20%, 30% and 50%). We start by 
using the line fits from Figure 3 to extract the redshift 
at which each sample of central galaxies crosses the se- 
lected value of f q . We then estimate the total stellar 
mass of central galaxies at this redshift using the derived 
relation of each number density sample (e.g., Equation 
2). Finally, we estimate a proxy for the group host halo 
mass using the total stellar mass and number density un- 
der Press-Schechter formalism (Press & Schechter 1974). 

Figure 4 shows the halo mass estimate of each number 
density sample at each value of f q as a function of red- 
shift (red points and lines). Also shown are predictions 
from two theoretical investigations of the halo quench- 
ing model. The green line represents the estimates from 
Dekel & Birnboim (2006), who follow an analytic ap- 
proach to study the formation of strong shocks in spher- 
ically symmetric halos. In this analysis, the critical mass 
value was estimated at log(M cr it/M©) ~ 11.7 with sig- 
nificant uncertainty of roughly half a dex. The magenta 
line is based on estimates from Cattaneo et al. (2006), 
who fine-tune the value of M cr i t in their semi-analytic 
simulations to best fit the color-magnitude distribution 
of 2 ~ 0 galaxies and the Lyman-break galaxy luminosity 
function at z ^ 3. These analyses use different tools to 
study the same model, according to which galaxies that 
reside inside a host halo more massive than some critical 
mass are efficiently quenched at low redshift (horizontal 
green and purple lines in Figure 4). At high redshift, 
high infall rates of cold narrow streams allow star forma- 
tion in higher mass halos (diagonal lines). The upturn 
redshift in Dekel & Birnboim (2006) is allowed in the 
range 1 < z crit < 3. 

Theoretical investigations of the halo quenching model 
agree with the observed result in that there seems to 
exist a threshold host halo mass above which star for- 
mation in central galaxies is suppressed efficiently. The 
implied typical range of the threshold halo masses 10 12 < 
M cr it/M© < 10 13 is consistent with the values favored by 
Cattaneo et al. (2006) and by Dekel & Birnboim (2006). 
We note that in addition to halo mass, local properties 
of the stellar component of galaxies (as traced by, e.g., 
central mass density and velocity dispersion) likely help 
determine star formation activity in these galaxies. 

4. ENVIRONMENTAL QUENCHING OF 
SATELLITE GALAXIES 

Measurements of the properties of individual satellites 
require that galaxies are uniquely associated to their par- 
ent halo, often relying on extensive spectroscopic data 
sets. Here we follow an alternative approach and perform 
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Figure 5. Demonstration of the technique used to derive quies- 
cent satellite galaxy fractions. Rest frame colors of satellite galaxy 
candidates were measured in apertures of size r v i r around selected 
centrals (solid red circle) and used to separate candidates into qui- 
escent and star forming bins. A similar measurement was per- 
formed around randomly selected positions in the imaging field 
(dashed red circle). The total number of quiescent sources in all 

background apertures ( N q 9 ) was subtracted from the total number 
of quiescent satellite candidates (7V| C ), and the result was divided 
by the difference between all background and satellite candidate 
galaxies. The resulting fraction ( fq at ) represents the average frac- 
tion of quiescent satellite galaxies around all centrals at a given 
redshift. While satellite galaxies cannot be identified individually 
around any one central, the ensemble of 6000 groups across all 
redshifts allows a clear statistical separation between satellite and 
background galaxies. 

a statistical analysis of the fraction of quiescent satel- 
lites around identified central galaxies. Using photomet- 
ric redshift measurements we quantify and subtract the 
contribution of foreground and background sources and 
derive the evolving quiescent satellite fraction in halos 
binned by mass and redshift. For the bulk of this analy- 
sis we utilize the n = 4x 10 -4 Mpc -3 sample in order to 
maximize completeness, redshift and source density. At 
this value of n, galaxies are at least 10 times more mas- 
sive than the 90% completeness limit at z < 1.8, allowing 
us to study the satellites of a mass evolving sample of 
central galaxies with a satellite-central mass ratio >1:10. 
At higher values of n completeness reaches 90% at lower 
redshift while at lower n the density of satellite galaxies 
is insufficiently high for background subtraction. 

4.1. Average Quiescent Fractions 

We select all satellite galaxy candidates from the Ul- 
traVISTA and 3D-HST/CANDELS catalogs within the 
virial radius of each central galaxy and within dz < 0.05 
and Msateiiite/Mcentrai > 0.1. Galaxies in SDSS were 
selected from the NYU-VAG catalog within dz < 0.01 
of their central galaxies. Virial radius estimates at all 
redshifts were acquired from the Guo et al. (2011) semi- 
analytic model. Redshift thresholds were defined based 
on the redshift accuracy of the catalogs and mass lim- 
its were determined to maximize the number of detected 
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Figure 6. Environmental quenching of satellite galaxies. Red 
points and blue diamonds represent the fraction of quiescent central 
and satellite galaxies, respectively. Black squares show the fraction 
of quiescent galaxies which are at the same redshift and mass as 
the satellites but are selected from randomly positioned apertures 
(“background” sources). The blue and black lines are fits to the 
data and the red line is the same as in Figure 3. Pale points 
represent the same measurements at stellar masses below the 90% 
completeness limit. Although all galaxies get quenched over time, 
central, satellite and background galaxies do so at different times 
and rates. While central galaxies are quenched early, star formation 
suppression in their satellites is delayed by a few Gyr. In addition, 
the quenching rate of satellites is faster than that of background 
galaxies at the same mass. The observed quenching rate and onset 
of satellite galaxies imply that star formation in these galaxies is 
strongly affected by environmental processes. 


satellites according to the calculated completeness value 
in stellar mass. 

We measure the rest frame UV J color (■ ugJ in SDSS) 
of each candidate satellite and count the number of quies- 
cent and star forming galaxies. We then perform an iden- 
tical measurement in ten randomly positioned apertures 
per galaxy to get an estimate of background and fore- 
ground source contamination. Finally, we subtract the 
average number of contaminating sources in each selec- 
tion category (quiescent or star forming) from the num- 
ber of candidate satellites in the same category. From the 
resulting distribution of UV J and ugJ values we calcu- 
late the average fraction of quiescent satellite galaxies at 
each redshift. Figure 5 shows a demonstration of this 
procedure. We also measured the fraction of quiescent 
background and foreground galaxies and the standard 
deviation of values in the ten stacks of randomly posi- 
tioned apertures. 

Figure 6 shows the fraction of quiescent galaxies as a 
function of lookback time and redshift at n = 4 x 10 -4 
Mpc -3 . Red data points represent central galaxies, black 
points show galaxies in randomly positioned apertures 
and blue points show the fraction of quiescent satellite 
galaxies. The pale pink data point at 2 = 3 and the gray 
and light blue data points at z > 2 reflect the same mea- 
surement where the 3D-HST/CANDELS catalog is less 


than 90% complete in stellar mass. The plotted blue er- 
ror bars represent the statistical measurement error and 
are equal to the standard deviation of measured values 
from the 10 random apertures in each redshift bin. The 
red error bars are the propagated \f~N noise estimates of 
the central quiescent fraction measurement. The black 
and blue solid lines show a linear fit to the data at 2 < 2 
and the red line shows the line fit to the central galaxy 
points in redshift-log(/ g ) space from Figure 3. 

4.2. Star Formation Quenching in Satellite Galaxies 

The most striking result that is evident from Figure 6 is 
the increase in the fraction of quiescent satellite galaxies 
compared to the overall background population. While 
on average all galaxies become more quiescent with time, 
they do so at different epochs. Central galaxies, which 
are selected to be more massive than their satellites, in- 
crease their quiescent fraction at early times. Satellites 
and background galaxies exhibit a mostly constant rate 
of increase in their quiescent fraction since z ~ 2. At that 
redshift, it is not possible to separate satellites from the 
general population of similar mass galaxies solely from 
their star formation activity. As time progresses, the qui- 
escent fraction of satellite galaxies increases compared to 
background galaxies and by z ~ 0, the fraction of non 
star forming satellite galaxies is nearly as high as that of 
their centrals. 

Satellite and background galaxies in our sample were 
selected to be identical in stellar mass and redshift, im- 
plying that the different quiescent fractions at z < 2 
are likely attributable to the environment in which these 
galaxies reside. In addition to halo quenching, which 
is expected to influence all galaxies, satellite galaxies in 
groups are subject to processes such as gas and stellar 
stripping, harassment and strangulation, all of which 
may assist in suppressing star formation. Moreover, 
satellite orbits in the host halo are typically eccentric 
(e.g., van den Bosch 1999) and as a consequence they 
spend a large fraction of the time outside of the virial ra- 
dius. This implies that the quenching time scale of each 
satellite galaxy depends on its orbit, as well as on any 
interactions that the galaxy experiences. Therefore, en- 
vironmental processes may act rapidly to suppress star 
formation in some satellite galaxies or may take a long 
time to quench other galaxies. Nevertheless, the aver- 
age quiescent fraction of the ensemble of satellite galax- 
ies increases over time compared to that of background 
galaxies at the same mass. 

4.3. Delayed Onset of Satellite Quenching 

Quenching of satellite galaxies in groups does not seem 
to begin immediately upon their infall into the group 
halo. Wetzel et al. (2012) utilized galaxy group cata- 
logs based on SDSS data in addition to a high resolution 
7V-body simulation in order to estimate the time delay 
between satellite galaxy infall and quenching. They find 
that star formation in satellites continues unaffected for 
2-4 Gyr after first infall. When quenching begins, it takes 
place on a rapid time scale that depends on the stellar 
mass of the satellites. Wetzel et al. (2012) conclude that 
satellite galaxy quenching is similar to that of central 
galaxies for 2-4 Gyr after being accreted into a group 
halo, and is significantly more rapid thereafter. 

Our results indeed show that satellite quenching is de- 
layed in respect to their centrals, although we cannot 
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Rapid central galaxy quenching 



Red shift ► 


Figure 7 . Top panel: star formation in central galaxies is sup- 
pressed (at least in part) by halo quenching, which acts after the 
halo reaches a critical mass of a few times 10 12 M©. In halos more 
massive than the critical mass, strong shocks are created at the 
group virial radius and inhibit subsequent star formation by heat- 
ing infalling gas. Mass and size growth continue in quenched galax- 
ies through mergers. Bottom panel: satellite galaxies in the same 
groups experience environmental quenching which starts a few Gyr 
after their centrals quench (see Figure 6). Environmental processes, 
such as gas and stellar stripping, harassment and strangulation, as- 
sist in shutting off star formation in satellite galaxies. While the 
quenching time scale of any individual satellite depends on its or- 
bit and on the interactions it experiences, the average quiescent 
fraction of all satellites increases with time compared to that of 
background galaxies at the same mass. 

measure the delay time scale as the current data hold no 
information regarding satellite accretion times. Instead, 
we estimate the initial quenching time delay by compar- 
ing the time difference between the onset of central and 
satellite galaxy quenching. To do so we measured the 
time difference between the central and satellite galaxy 
curves at f q ~ 0.1 (horizontal distance between the red 
and blue curves in Figure 6), when the effects of environ- 
mental processes first separate the quiescent fraction of 
satellites from that of background galaxies. We repeated 
the measurement for an additional number density selec- 
tion, n = 6x 10 -4 Mpc -3 , by finding the best fit curves to 
the data in the same way as discussed in Section 4.1 and 
presented in Figure 6. We note that selected satellites 
at this value of n reach M*/M© = 6.5 x 10 9 at z = 0, 
roughly twice the stellar mass of the Large Magellanic 
Cloud. 

The derived initial time delay is 1.5 ±0.3 Gyr and 
2.2±0.5 Gyr for the n = 4 x 10 -4 Mpc -3 and n = 6 x 10 -4 
Mpc -3 samples, respectively. Both values are slightly 
shorter than the estimate of 2-4 Gyr that was calculated 
for similarly massive galaxies at z = 0 from Wetzel et al. 
( 2012 ). 

4.4. Sources of Uncertainty 

The results that are presented in this study are sub- 
ject to potential biases which may have been introduced 
by our galaxy selection approach. For example, we se- 
lect satellite galaxies within a fixed mass threshold from 


their group centrals. Since we do not directly measure 
how the mass of satellite galaxies evolves over the stud- 
ied redshift range, we cannot predict the corresponding 
actual evolution in the mass ratio value. An alterna- 
tive selection which includes all satellite galaxies above 
a fixed threshold mass, would imply an increasing mass 
ratio with time and a large number of low mass satellites 
at low redshift. Consequently, the contribution of lower 
mass satellites may affect the measured quiescent frac- 
tions of satellite and background galaxies alike and may 
change the rate of their evolution over time. We discuss 
this selection bias further in Appendix C. 

An additional assumption that may affect our mea- 
surements if sufficiently inaccurate is that of a constant 
population of satellite galaxies in the studied halos. If 
a large fraction of satellites at any given redshift are 
subject to accretion and merging, our measurement of 
satellite quenching may not be dominated by environ- 
mental processes as satellites do not spend enough time 
in the vicinity of the group for such processes to take 
effect. However, we note that recent studies find that 
most satellite galaxies orbit the halo on eccentric orbits 
and that they spend a large fraction of their time out- 
side of the group halo virial radius (but are still affected 
by the group; e.g., Tal et al. 2013; Wetzel et al. 2013). 
This implies that even in the presence of accretion and 
merging, the overall population of satellite galaxies may 
be mostly stable over time. 

Finally, the selection of constant number density of 
central galaxies may have affected the measured quies- 
cent galaxy fractions. Even under the assumption that 
overall, a fixed cumulative number density selection truly 
follows a given population of galaxies over time, this tech- 
nique may inadequately follow populations of only quies- 
cent or star forming galaxies. If instead we effectively fol- 
low vastly different selections of galaxy populations, our 
interpretation may suffer from significant biases. Never- 
theless, since cumulative number density is probably our 
current best approach for tracing galaxies over different 
epochs, we utilize it throughout this study. 

5. SUMMARY 

Galaxy environments have long been postulated to 
play a crucial role in regulating the star formation rates 
of their resident galaxies. Here we studied the link be- 
tween environment and star formation in groups by an- 
alyzing the fraction of quiescent central and satellite 
galaxies in the redshift range 0 < z < 2.5. We followed 
central galaxy populations at fixed values of cumulative 
number density and traced group satellites with mass ra- 
tios as low as 1:10. In addition, we analyzed background 
galaxies (which act as centrals in their own host halos) 
at the same mass and redshift as the satellites. Cen- 
tral galaxies in our sample span a wide range of stellar 
masses (3.5 x 10 9 < M/M© < 2.2 x 10 11 ) over seven 
redshift bins and four number density selection values. 
We found clear evidence for star formation suppression 
in group halos and measured different quenching times 
for satellite and central galaxies. 

Initial star formation suppression in central galaxies 
takes place early. Figure 3 shows that at z ~ 2.5, at 
least 30% of the centrals in each of the two most massive 
samples are already quenched. Even the lowest mass 
centrals in our sample had already started quenching by 
z ~ 1.5. We utilized these measurements to calculate a 
threshold halo mass for star formation suppression and 
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Central galaxy selection criteria for all number density samples 
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Selection 

Data set 

Redshift 

Central mass range 

A7 med 

(M e ) 

Rvir 

(kpc) 

Ncen 


SDSS 

0.02 < 2 : < 0.04 

10.70 <log {M/Mq) < 11.00 

6.2X10 10 

220 

667 

n — 6 x 10 -4 Mpc -3 


0.2 < 2 : < 0.4 

10.54 <log (M/M 0 ) < 10.85 

4.9xl0 10 

200 

301 


UltraVISTA 

0.4 < 2 : < 0.6 

10.45 <log(M/M 0 ) < 10.75 

4.0xl0 10 

180 

531 



0.6 < z < 0.85 

10.35 <log(M/M 0 ) < 10.65 

3.1X10 10 

150 

1251 

log (M*/M©) = —0.07 z 2 - 0.37z + 10.81 


0.85 < z < 1.2 

10.20 <log (M/M©) < 10.50 

2.2X10 10 

110 

2858 



1.2 < z < 1.8 

9.93 <log (M/M©) < 10.23 

1.2X10 10 

85 

970 


3D-HST 

1.8 < z < 2.5 

9.55 <log {M/Mq) < 9.85 

5.0X10 9 

60 

1827 



2.5 < 2 : < 3.5 

8.87 <log(M/M 0 ) < 9.17 

1.0x10® 

30 

3530 


SDSS 

0.02 < ^ < 0.04 

10.83 <log (M/Mq) < 11.13 

7.9xl0 10 

250 

251 

n — 4 x 10 -4 Mpc -3 


0.2 < 2 : < 0.4 

10.69 <log (M/Mq) < 10.99 

6.9xl0 10 

220 

261 


UltraVISTA 

0.4 < z < 0.6 

10.64 <log (M/Mq) < 10.94 

6.1xl0 10 

200 

460 



0.6 < z < 0.85 

10.57 <log (M/Mq) < 10.87 

5.2xl0 10 

180 

1080 

log (M*/M®) = — 0.10z 2 - 0.18z + 10.90 


0.85 < z < 1.2 

10.45 <log (M/Mq) < 10.75 

4.1xl0 10 

150 

2281 



1.2 < 2? < 1.8 

10.25 <log (M/Mq) < 10.55 

2.5xl0 10 

110 

707 


3D-HST 

1.8 < 2? < 2.5 

9.92 <log(M/M 0 ) < 10.22 

1.2xl0 10 

75 

915 



2.5 < z < 3.5 

9.29 <log(M/M©) < 9.59 

2.8xl0 9 

40 

3054 


SDSS 

0.02 < 2 : < 0.04 

10.97 <log (M/M q ) < 11.27 

l.lxlO 11 

350 

50 

n = 2 x 10 — 4 Mpc -3 


0.2 < 2? < 0.4 

10.89 <log (M/M q ) < 11.19 

l.OxlO 11 

300 

165 


UltraVISTA 

0.4 < 2? < 0.6 

10.87 <log (M/Mq) < 11.17 

9.8xl0 10 

270 

250 



0.6 < 2 : < 0.85 

10.83 <log (M/Mq) < 11.13 

9.0xl0 10 

240 

610 

log(M*/M®) = — O.IOz 2 - 0.04z + 11.03 


0.85 <z<1.2 

10.76 <log(M/M 0 ) < 11.06 

7.7xl0 10 

200 

1294 



1.2 < 2 : < 1.8 

10.61 <log (M/Mq) < 10.91 

5.6xl0 10 

150 

526 


3D-HST 

1.8 < z < 2.5 

10.35 <log (M/Mq) < 10.65 

3.2X10 10 

100 

595 



2.5 < 2 : < 3.5 

9.86 <log (M/Mq) < 10.16 

9.9x10® 

60 

1601 


SDSS 

0.02 < 2 : < 0.04 

11.05 <log (M/M q ) <11.35 

1.5X10 11 

380 

21 

n = lx 10 -4 Mpc -3 


0.2 < 2 ? < 0.4 

11.02 <log (M/M q ) <11.32 

1.4X10 11 

350 

133 


UltraVISTA 

0.4 < 2 : < 0.6 

11.00 <log (M/Mq) < 11.30 

1.4X10 11 

330 

167 



0.6 < 2 : < 0.85 

10.99 <log(M/M 0 ) < 11.29 

1.3X10 11 

310 

369 

log (M*/M®) = — 0.08z 2 - 0.02z + 11.17 


0.85 < z < 1.2 

10.96 <log(M/M e ) < 11.26 

1.2X10 11 

250 

714 



1.2 < z < 1.8 

10.85 <log (M/M q ) < 11.15 

9.2X10 10 

190 

256 


3D-HST 

1.8 < z < 2.5 

10.65 <log (M/Mq) < 10.95 

6-lxlO 10 

140 

437 



2.5 < 2 : < 3.5 

10.27 <log(M/M 0 ) < 10.57 

2.6xl0 10 

70 

759 


found it to be in the range 10 12 < M crit /M® < 10 13 
(Figure 6). This value is relatively constant at 1.3 < 
z < 3.0 and it is consistent with theoretical predictions 
from the halo quenching model. Local processes, which 
were not considered in this analysis, likely contribute to 
quenching as well. 

Satellite galaxies begin quenching a few Gyr after the 
onset of star formation cessation in the centrals. Figure 
6 shows that at z ~ 1.5 the average fraction of quies- 
cent satellite and background galaxies is similar to one 
another but significantly lower than that of centrals. At 
that redshift, centrals had already been quenching for 
1-3 Gyr. At z ~ 0, the fraction of quiescent satellites 
approaches and nearly reaches the fraction of quiescent 
centrals. This delayed quenching was also reported by 
Wetzel et al. (2012) in low redshift groups. 

Additional evidence for the difference between central 
and satellite galaxy quenching comes from the different 
quiescent fractions of satellite and background galaxies at 
z < 1.5. Since sources that were selected as background 
galaxies act as centrals in their respective halos, the black 
and blue lines in Figure 6 represent the difference be- 


tween central and satellite quenching at the same mass 
and redshift. This is in agreement with studies who find 
that satellite galaxies are on average more quenched than 
same-mass centrals (e.g., Wetzel et al. 2012). Overall, 
the observed quenching in satellites is consistent with en- 
vironmental processes (mainly gas heating and stripping) 
which act to suppress star formation in galaxy group ha- 
los. 

In summary, we find that while star formation gets 
suppressed in all galaxies in groups, the processes that 
govern this suppression are not universal (Figure 7). Our 
analysis shows that central galaxies in groups undergo 
quenching which is consistent with the halo quenching 
model while their satellites experience additional quench- 
ing by environmental processes. The onset of star forma- 
tion suppression also varies between centrals and satel- 
lites with early central quenching and delayed quenching 
of satellites. 


We thank Guillermo Barro and Alis Deason for engag- 
ing discussions and for commenting on an early draft of 
the paper. 
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Figure 8. Rest frame color based selection of quiescent galaxies from SDSS, UltraVISTA and 3D-HST/CANDELS. Gray points represent 
the distribution of all galaxies in the catalogs while red and blue points show selected central galaxies and satellite candidates, respectively. 
At z > 0.2 quiescent galaxies were selected from a rest frame UV J diagram according to the thresholds found by Muzzin et al. (2013a). 
Lower redshift galaxies were selected in rest frame ugJ space to roughly match the high redshift selection. The enhanced vertical scatter 
in the SDSS points largely results from relatively poor u — band photometry. 
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APPENDIX A 

CUMULATIVE FIXED NUMBER DENSITY 
SELECTION PARAMETERS 

Central galaxy selection criteria are presented in Table 
2. The selection column shows the fixed cumulative num- 
ber density value of each sample, as well as its derived 
mass-redshift relation (Section 2.2). Stellar mass ranges 
were calculated such that each log (M/M©) bin of size 0.3 
dex had a median central galaxy mass according to its 
respective mass-redshift relation. Virial radius estimates 
were calculated using catalogs based on the semi-analytic 
models of Guo et al. 2011 for galaxies at the same red- 
shift and mass. The number of central galaxies in each 
bin are given in the final column of Table 2. 

APPENDIX B 

UVJ AND ugJ SELECTION 

Quiescent galaxies at z > 0.2 were selected from the 
UltraVISTA and 3D-HST/CANDELS catalogs based on 
their location on a rest frame UV J diagram. The ancil- 
lary data sets and SED modeling techniques that were 
used to derive these colors for both catalogs are similar 
and therefore so are the galaxy selection limits that were 
utilized in this study. The middle and right-hand panels 
of Figure 8 show the distribution of galaxy colors in the 
catalogs, as well as the selection thresholds (black lines) 
which were adopted from Muzzin et al. (2013a). 

At z ~ 0, rest frame color estimates require very small 
corrections from observed colors. More specifically, u — g 
color values in the NYU-VAG catalog were derived rely- 
ing heavily on u — band photometry from SDSS. The re- 
sponse function of the SDSS photometric system in this 
filter band is notoriously poor and as a consequence flux 
estimates suffer from large errors (Fukugita et al. 1996). 
While the u — g and U — V colors that were used in this 
study are very similar in wavelength coverage, the en- 
hanced scatter in u — band photometry implies that the 
u — g measurement is more uncertain. The left panel of 
Figure 8 shows the distribution of galaxy colors in SDSS. 
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It is evident from the figure that the division between star 
forming and quiescent galaxies is far from perfect as the 
two populations are somewhat mixed. 

In order to match the selection in ugJ to that of UV J 
we followed Williams et al. (2009) and Muzzin et al. 
(2013a) and tried to contain as much of the quiescent 
galaxy peak while minimizing contamination from star 
forming galaxies. The lines in the left hand panel of Fig- 
ure 8 show these selection thresholds. Despite significant 
scatter, the majority of all quiescent central and satellite 
galaxies are likely included in our selection. 


APPENDIX C 

CONSTANT MASS SELECTION OF CENTRAL 
AND SATELLITE GALAXIES 

One potential bias that could affect our analysis of the 
evolving quenched fractions of satellite galaxies (Section 
4) stems from the number density selection approach that 
we utilized in this study. Since we do not directly follow 
the mass growth of satellite galaxies but rather require a 
fixed ratio between central and satellite masses, some of 
the implied evolution of the quiescent fractions may be 
the result of mass growth over the same redshift range. 
Here we test the robustness of our results by repeating 
the analysis for a constant mass selected sample of central 
and satellite galaxies. Central galaxies at all redshift s are 
in the mass range 10.75 < log (M/M q) < 11.05 while all 
satellite and background galaxies are more massive than 
log(M/M 0 ) > 9.75. Figure 9 demonstrates that even 
in the absence of any mass evolution, satellite galaxies 
experience delayed quenching after their centrals, at a 
more rapid rate than background galaxies at the same 
mass and redshift. We conclude that the observed trends 
that are presented in this study do not entirely result 
from the mass evolution of satellite galaxies. 
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